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We propose that the concept of multidimensional residues can be used to directly extracting the 
coefficients of scalar master integrals (with single propagators only) from one-loop Feynman integrals 
with generic power of propagators. Unlike the usual integration-by-parts (IBP) technique, where 
one has to solve iteratively a complicated set of equations to carry out the reduction and determine 
the coefficients of scalar master integrals, using multidimensional residues provides the possibility 
of directly extracting the coefficients of the master integrals. As the first application of this idea, 
we show how to directly extract the scalar box integral coefficients. 
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Introduction. Scattering amplitudes are one of the 
most widely studied quantities within the framework of 
perturbative quantum field theory. They are important 
ingredients for computing various observables that can 
be used to test the theory in which they are defined. 
Studying scattering amplitudes can also shed light on 
potential hidden symmetries of the defining theory that 
are not clearly visible from the Lagrangian. 

From both phenomenological and theoretical perspec- 
tives, one needs to study scattering amplitudes beyond 
tree-level. In general, such amplitudes are written in 
terms of tensor Feynman integrals, which can then be 
reduced down to a set of irreducible integrals. At one- 
loop level, these irreducible integrals comprise scalar in- 
tegrals with no more than four propagators. This fact 
has played an important role in recent developments in 
computing one-loop scattering amplitudes using the so- 
called on-shell approach [ll-Q (for recent reviews on this 
approach, see e.g. [8, §]). 

The on-shell approach provides an efficient means of 
constructing the one-loop amplitude from its singularity 
structures. One can decompose the one-loop amplitude 
into: the cut-constructible part, i.e. the scalar four-, 
three-, two- and one-point scalar integral functions multi- 
plied by rational coefficients; and the rational terms. The 
scalar integral functions contain all (poly) logarithmic de- 
pendence of the amplitudes. Hence by computing and 
matching the branch cuts of the amplitude and of its 
decomposition in terms of the scalar integral functions, 
one can extract information on the rational coefficients of 
the scalar integral functions. For instance, the coefficient 
of a scalar box integral can be determined by imposing 
a quadruple cut that uniquely isolates the contribution 
of this scalar box integral j3(]. A quadruple cut condi- 
tion completely freezes the loop momentum (if it is taken 
as four-dimensional), therefore matching the cuts deter- 
mines the coefficient of the scalar box integral in a com- 
plete algebraic way. Extracting the coefficients of scalar 
three-, two- and one-point integrals is more complicated, 
as in these cases cutting the propagators of the desired 
integral topology does not uniquely isolate the desired 
coefficient. To determine the desired integral coefficient. 



one has to disentangle the contribution of other integrals 
involving not only the cut propagators but also extra 
propagators. This can be done following e.g. Refs. 0,0). 
Since the scalar integral functions are well-known, after 
determining all coefficients of these integral functions, 
one has the result for the cut-constructible part of the 
amplitude. The second part, i.e. the rational terms can 
not be computed from four-dimensional cuts, since they 
are purely rational functions and do not have any cuts. 
However, when moving to /^-dimensions, these rational 
terms also develop branch cuts and can thus be computed 
from the Z)-dimensional cuts 0, • Another possibility 
of computing the rational terms is to use on-shell re- 
cursion relations inspired by the BCFW recursions [lH 
for tree- level amplitudes Combining the cut- 

constructible part and the rational terms, one arrives at 
a form of the one-loop amplitude in terms of scalar master 
integrals without performing an explicit tensor reduction. 

In Ref. the reduction of one-loop amplitudes was 
carried out at the integrand level. By investigating the 
structure of the integrand of a generic one-loop ampli- 
tude, the authors of Ref. [sj were able to decompose the 
integrand of a one-loop amplitude into integrands cor- 
responding to scalar master integrals and the spurious 
terms, where the spurious terms depend on the loop mo- 
mentum and yield vanishing contribution upon integra- 
tion. The coefficients of scalar master integrals as well as 
of the spurious terms can then be extracted by imposing 
four, three, two and one on-shell conditions on the inte- 
grand of the one-loop amplitude and its decomposition, 
and then matching them. The rational terms of the one- 
loop amplitude can also be computed using the method 
proposed in Ref. 

The above procedure of directly extracting the scalar 
master integral coefficients applies to generic tensor one- 
loop Feynman integrals with single propagators. If there 
is higher power of propagators involved, one has to resort 
to other techniques. One of the most common such tech- 
niques is to use the so-called integration-by-parts (IBP) 
identities The IBP technique allows one to derive 

a set of equations relating Feynman integrals with high 
powers of propagators to those with lower powers. Solv- 
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ing these equations, one is able to reduce Feynman inte- 
grals with high powers of propagators. This reduction 
procedure is an iterative one. By deriving and solv- 
ing the IBP equations iteratively, one finally achieves a 
complete reduction of one-loop Feynman integrals with 
generic power of propagators down to master integrals 
with single propagators only. The procedure of solving 
the IBP equations can be rather complicated, and effi- 
cient algorithms are usually needed. 

In this paper, we propose that the concept of multidi- 
mensional residues provides a new efficient means of di- 
rectly extracting the coefficients of scalar master integrals 
from one-loop Feynman integrals with generic power of 
propagators. In contrast to the IBP method, using multi- 
dimensional residues allows one to directly determine the 
coefficients of scalar master integrals without iteratively 
solving the IBP equations. Multidimensional residues 
have played an important role in complex analysis and 
geometry, but are so far not widely used in physics. There 
have been some physics references discussing their appli- 
cations in string theory [l5l | and in studying leading sin- 
gularities of scattering amplitudes in Af = A super Yang- 
Mills theory 16|. We argue that when generic power 



of propagators is involved in a Feynman integral, mul- 
tidimensional residues have the potential to efficiently 
extract the coefficients of master integrals. As the first 
application and test, we show how to directly extract the 
coefficients of scalar box integrals appearing in the reduc- 
tion of one-loop Feynman integrals with generic power of 
propagators, using multidimensional residues. 

Before we proceed, we first give a brief introduction to 
the concept of multidimensional residues. 

Multidimensional residues. The discussion in this 
section is essentially following Ref. [l3|- 

Suppose one has n holomorphic functions gi, 32, • " ' 5n 
defined in an open set U C C", where the n functions 
have a single common zero p in U, one can associate to 
any holomorphic function h G 0{U) the local residue at 
p of the memomorphic form 



/i(x) 



5i(x)52(x)---g„(x) 

This then defines a C-linear operator 

h{x.)dx 



with dx. = dxif\dx2/\- ■ ■ /\dXn- 

(1) 
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5i(x)---5n(x)/ ' (27ri)" 7rg(£)' 

(2) 

where rg(e) is the real n-dimensional cycle 

rg(e) = {xe[/: |5,(x)| -e„^ = l,...,n} (3) 

with orientation defined by the n-form darg(gi) A • ■ • A 
darg((7„), and e — (ei,...,e„) is any n-tuple of suffi- 
ciently small, positive real numbers. 

As a straightforward generalization of the Cauchy for- 
mula for one complex variable, the residue of the mem- 



omorphic form uj with 
given by 

/i(x)dx 



[x. 
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[xi -PiY 



• • • a„! \dx\^ ■ ■ ■ dxn" ' 



ai 



(4) 

This relation will be particularly useful when we discuss 
below the extraction of the coefficients of scalar master 
integrals arising in the reduction of one-loop Feynman 
integrals with generic power of propagators. 

Scalar one-loop integral coefficients from multi- 
dimensional residues. To illustrate how the idea of 
multidimensional residues can be applied to determining 
the scalar one-loop integral coefficients, let us start from 
one-loop amplitudes with single propagators only. Ac- 
cording to Ref. [5] , the integrand of any m-point one- loop 
amplitude can be decomposed as 



N{q) 



E 



d{ijkl) + d(q; ijkl) 



DoDi ■ ■ ■ Dra-i D.DjDkDi 

i<]<k<l ■' 

E™"^ c{ijk) c{q] ijk) h{ij) + ij) 



i<j<k 



i<j 



E 



a{i) + d{q; i) 



(5) 



where Di — {q+PiY —mf; d{ijkl),c{ijk),b{ij),a{i) are 
the coefficients of the scalar 4-, 3-, 2-, 1-point integral, re- 
spectively; d{q; ijkl), c{q; ijk), b{q; ij), d{q; i) are the spu- 
rious terms. We do not write down the spurious terms 
without denominators (propagators) in the above expan- 
sion, that is, we will always work in a renormalizable 
theory and choose the renormalizable gauge where such 
terms vanish. 

As far as the coefficients of scalar one-loop integrals 
are concerned, we can take the loop momentum as four- 
dimensional. We choose the same decomposition of the 
loop momentum as in Ref. Q 



-Pi 



(6) 



where are four linearly independent massless momenta 
constructed from the external momenta. We will view Xi 
as independent complex variables. The integrand of the 
amplitude then becomes a function of these complex vari- 
ables. The scalar box integral coefficient can be deter- 
mined by computing and matching the residues of both 
sides of Eq. ^ at the pole determined by the on-shell 
conditions 



Do = Di = D2 = D3 = 0. 



(7) 



It is convenient to make a change of variables from Xi 
to Di when computing the residue at Dq = Di = D2 = 
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Da = 0. This change of variables introduces a Jacobian 



Ja = dct 



/ d{D„,Di,D2,D3) \ 
\ 9(a;i,X2,a;3,a;4) / 



(8) 



into the residue. However, for an amphtude with sin- 
gle propagators, this Jacobian cancels out when match- 
ing the residues of both sides of Eq. (O . The remaining 
terms yield exactly the same relation between the numer- 
ator N{q), the scalar box coefficients and the 4-point like 
spurious terms as that obtained in Ref. Q . The relations 
to determine the scalar triangle, bubble and tadpole co- 
efficients can be obtained by viewing one, two and three 
of Xi as input parameters, then computing and matching 
the residues of the amplitude and its decomposition at 
the poles determined by the on-shell conditions isolating 
the desired scalar integral coefficients. 

We now argue that the idea of multidimensional 
residues can also be used to directly determining the co- 
efficients of scalar one-loop integrals arising in the reduc- 
tion of one-loop Feynman integrals with generic power of 
propagators. The main point of this argument is that at 
the poles of the propagators, the integrand of the origi- 
nal integrals with high power of propagators has the same 
(multidimensional) residues as the integrand of their de- 
composed form in terms of scalar master integrals with 
single propagators. As the first application and test of 
the idea, we show how to directly determine the scalar 
box integral coefficients using multidimensional residues. 
As in the reduction of Feynman integrals with single 
propagators, there exist contributions of spurious terms 
at the integrand level. The 4-point like spurious terms 
now look like 



d{q; ijkl, hjhjhkhi) 



(9) 



where the numerator has the same structure as d{q; ijkl) 
in Eq. ([SJ , but now we can have more power of propaga- 
tors. 

Following the idea proposed above, the relation to de- 
termine the scalar box coefficient from multidimensional 
residues is as following: 



ReSr, 



N{q)d^x 



r)9o nsi . . . r)yn.-i 
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d(ijkl) 
D,D,DkD 



E 



i<j<k<l 
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(10) 



where the subscript p denotes the pole determined by 
the on-shell conditions Di = Dj — Dk = Di = 0. In 
the following we take i = 0,j = l,k = 2,1 = 3 as an 
illustrative example, the procedure for other cases are 
completely analogous. 



As mentioned above, the numerator of the 4-point like 
spurious terms has the same structure as in Ref. jBf, 
namely 

d{q; 0123, hohMa) = d{0123, hohih2h3)T{q), (11) 
where 

T{q)^TT[i^+P^)tj2hl5], (12) 

hjh are two of the four massless momenta constructed 
from the external momenta as 



^1 + aih 



ki 

'3 



k2—l2+ Ci2h, 



(13) 



with k^ = p^-pQ. 

We first write the the propagators Di{i — . . . 3) in 
terms of the complex variables Xj (j = 1 . . . 4) as 

-Do = lixiX2 - 4^x3X4) - do, 
Di - Do = do - di + -f{xiai + X2), 
D2 - Do = do - d2 + ^{x2a2 + a;i), 
D3^Do = do- d3 + 2[xi(k3 ■ h) + X2{k3 ■ h) 

+ X3ik3 ■ I3) + X4{k3 ■ U)i (14) 

where di — mf — kf . From these equations one can solve 
Xj in terms of Di and compute the Jacobian in Eq. . 
Using Eq. one obtains the residue on the l.h.s. of 
Eq. (fTO]) at the pole as 



.90! • • - 53! 



,h{q^{x,{D,))) 



dD^o" 



■dDf 



0,1,2,3—0 

(15) 



where gi — Qi ~ ^, qo correspond to the two solutions 
Eq. (HID. Ri{q^) is 



to 



^4(90^) 



n 



15^0,1,2,3 



(16) 



and J4 is the Jacobian introduced in Eq. ([8]). Before 
computing the residue on the r.h.s. of Eq. ([T0| . we 
point out that, by an explicit computation of the Ja- 
cobian J4, we find that it has the same Xi dependence 
as the 4-point like spurious terms. Therefore the ratio 
d{q; 0123, hohih2h3) / J4 becomes independent of Xi and 
thus independent of Di. Consequently all spurious terms 
with any of the propagator power higher than 1 yield 
vanishing contribution to the residue, since they would 
involve a derivative with respect to at least one of Di. 
The only remaining contributions to the residue are 



rf(0123) /3J(0123,1111) 
4^5 



(17) 
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where /? = 1/(1 — aia2),7 = 2Zi ■ Z2. Note that as the 4- 
point hke spurious term in the case of single propagators, 
Jiilo) = ^^4(90")- Matching the residues on both sides 
of Eq. ([TUl) then gives the scalar box integral coefficient 



d(0123) 

-Do, 1,2, 3—0 

(18) 

which obviously reproduces the scalar box integral coef- 
ficient for the case where all propagators are single prop- 
agators. We have checked with several Feynman inte- 
gral examples involving different power of propagators 
that, the scalar box integral coefficient computed from 
Eq. (|18l) exactly agrees with that obtained from the IBP 
technique [l^. It is reasonable to expect that, as in the 
case of single propagators, the above procedure of using 
multidimensional residues to directly extract the scalar 
integral coefficients continues to work for triangle, bub- 
ble and tadpoles, provided that one obtains all the co- 
efficients of spurious terms or finds a way to avoid solv- 
ing them. We will explore this in a forthcoming paper. 
Moreover, the idea presented here may apply to high loop 
level. We will also explore this point in the near future. 

Conclusions. In this paper we proposed a new efficient 
means of directly extracting the scalar one-loop integral 
coefficients from one-loop Feynman integrals with generic 
power of propagators, using multidimensional residues. 
In contrast to the IBP technique, which is widely used for 
the reduction of Feynman integral with repeated propa- 
gators, using multidimensional residues allows to directly 
extract the coefficients of the master integrals without 
solving iteratively a set of IBP equations. As a first ap- 
plication and test of the idea, we showed how to directly 
extract the scalar box integral coefficients. 
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